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Abstract. Long-slit optical emission-line spectra of the H2O megamaser galaxies ESO 1Q3-G035, TXS 2226-184, 
and IC 1481 are evaluated in order to look for characteristics typical for water-megamaser galaxies. We present 
rotation curves, line ratios, electron densities, temperatures, and H/3 luminosities. The successful line-profile 
decompositions rest on d-Lorentzians with an additional parameter d to adjust the wings, rather than Gaussians 
or Lorentzians as basic functions. No significant velocity gradient is found along the major ax;is in the innermost 
2kpc of TXS 2226-184. IC 1481 reveals a spectrum suggestive of a vigorous starburst in the central kiloparsec 
10* years ago. None of the three galaxies shows any hints for outfiows nor special features which could give clues 
to the presence of H2O megamaser emission. The galaxies are of normal Seyfert-2 (ESO 103-G035) or LINER 
(TXS 2226-184, IC 1481) type. 
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1. Introduction 

Extragalactic H2O megamaser sources are orders of mag- 
nitude more powerful than strong A=1.3cm Galactic 
^ ' masers of the same type and have been discovered in about 
k> three dozen galactic nuclei, almost all of them of Seyfert- 
^ 2 or LIN ER (low-ionization nuclear emission-l i ne reg ion) 
5^ type (e.g. lBraatz et al.l l)l996|) : lGreenhill et alJ l|2003|i b 

According to the so-called unifi ed model 
ijAntonucci fc Milleii Il985t lAntonucci Il993|) . Seyfert-2 
objects contain a nuclear molecular torus seen edge-on 
that surrounds an X-ray luminous, winds and jets 
producing active galactic nucleus (AGN). It is plausible 
that this geometry is favorable for megamaser activity, 
because the necessary large column densities of warm 
(T > 400 K) and dense (n(Il2) > 10''cm~^) molecular gas 
enriched with H2O molecules are likely to be supplied. 
Interaction of molecular gas with a radio jet may also 
cause c onditions favorab le for the occurrence of H2O 
masers (|Peck et al.l l2003l) . 



Send offprint requests to: Nicola Bennert 

* Based on observations collected at the European Southern 
Observatory, Chile 
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This paper is the second one in a series ana- 
lyzing the spectral properties of water-vapor mega- 
maser galaxies at optic al wavelengths. In the first one 
l|Schulz &: Henkeil2003tl . emission-line profiles of IC 2560, 
NGC1386, NGC1052, and Mrkl210 were evaluated. 
Galactic rotation and outflow of narrow-line gas are com- 
mon features of this sample of water-megamaser galaxies. 
All decomposed line-systems exhibit AGN typical line ra- 
tios. For NGC 1052 and Mrkl210, maser emission trig- 
gered by the optically detected outflow components can- 
not be ruled out. 

In this work, we continue to explore the structure 
of galactic nuclei containing H2O megamaser sources 
(ESO 103-G035, TXS 2226-184, IC 1481) in order to look 
for distinguishing characteristics and to obtain clues to 
the nuclear geometry. 

1.1. Individual objects 

[Marshall etlll l)l979() and IPiccinotti eHIl l)l982l) classi- 
fied ESO 103- G035 (IRAS 183 3 3-652 8) as a Seyfert-1.9 
galaxy whereas iMorris fc WardI l)l988f) did not flnd any 
hints of broad line wings. Observations with EXOSAT 
(European Space Agency's X-ray Obser vatory) revealed 
variable, strong absorption of soft X-rays l) War wick et all 
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ll988l) . lHeisler fc Vadeil lll994l) find ESO 103-G035 among 
the warmest far-infrared emitters in their sample of IRAS 
(Infrared Astronomical Satellite) galaxies with spectral 
energy distributions peaking near 60 microns. It unveals 
a high excitation emission-line spectrum su perposed on a 
red st ellar continuum with absorption lines. iBraatz et al.l 
l|l996() discovered II2O megamaser emission in this galaxy. 

So far, the most luminous known II2O maser is 
found in the galaxy TX S 2226-184 (IRAS F22265-1826; 
iKoekemoer et all l|l995j) V Hubble Space Tele scope (HST) 
and V ery Large Array observations from iFalcke et alJ 
l(200r^ suggest that it is a galaxy of type later than SO with 
an inclination of ~70°. Their Ha + [Nil] map exhibits 
a gaseous, jetlike structure perpendicular to the galaxy's 
major axis and its nuclear dust lane. The 8.4 GHz radio 
continuum map shows emission that is elon gated in the 
same direction as the Ha + [N il] emission. iFalcke et al.l 
1I2OOOI) concluded that the nuclear accretion disk, the ob- 



scuring torus, an d the large-scale mo lecular gas layer are 
roughly coplanar. iTavlor et all l|2002l) found H I in absorp- 
tion towards TXS 2226-184, consisting of two components 
with widths of 125 km s""'^ and 420kms~^. They suggest 
that the H l absorption is produced in the central parsecs 
of the galaxy, on a scale similar to that which gives rise to 
the water maser emission. 

Little is known about IC 148 1 (IRAS 23168+0537 ) 
which was cl a ssified as LINER bv iHuchra eTall l)l982|) . 
iBraatz et al.l l(l996l) discovered its water megamaser 

emission. D ue to its amorphous appearance (e.g. 

Falcke et al.l (l200lf)') i t s incli nation is hard to determine: 



De Vaucouleurs et al ] lll99lh (RC3) list an inclination of 



i ^ 70° whereas i ^ 3 0° is given in IBraatz et al ] 11997). 
IVan den Bersh et all l|2002t) classify IC 1481 in their cat- 
alogue of host galaxies of supernovae as a peculiar Sb 
galaxy. 



2. Observations and data reduction 

The spectra described here were obtained by H. 
Domgorgen in May 1996 using the Boiler & Chivens 
spectrograph attached to the Cassegrain focus of the 
European Southern Observatory (ESO) 1.52 m telescope. 
Observations were made in two spectral ranges (3400- 
5400 A and 5400-7400 A) through the nucleus along the 
major axis of each galaxy with exposure times ranging 
between 1800 s and 3600 s. A log of the observations is 
given in Tabled The detector used was La Silla CCD No. 
39 (Loral Lesser) with 15 /im wide square pixels. The spa- 
tial resolution element is 0'.'68pix~^. Seeing and telescope 
properties limited the spatial resolution to the range 1'.'2 — 
2'.'5 which was determined by the width of spectra of stan- 
dard stars on the focal exposures. The 2" wide slit projects 
to a spectral resolution of ^ 2.7 A 130kms~^) as is con- 
firmed by the full width at half-maximum (FWHM) of 
comparison lines and the [O i] A6300 night-sky line. 

Standard reduction including bias subtraction and 
fiat-field correction was performed using the ESO 



MIDAS^ software (version Nov. 98). Night-sky spectra 
"above" and "below" any notable galaxy emission were 
interpolated in the region of the galactic spectrum and 
subtracted in each case. Excellent flux interpolation was 
achieved by rebinning the spectra to a scale of 0.97 A per 
pixel during wavelength calibration. The curve of iTiia 
II1977D was used to correct for atmospheric extinction. 
The spectra were flux calibrated using the standard stars 
LTT7987 or CD-32°9927. 

From the frames cleaned in this way, single rows were 
extracted (each row corresponding to 0'.'68 and 2" along 
and perpendicular to the slit direction, respectively). 
Along the "spatial axis" of the CCD, we identified the 
"photometric center" (that we choose as the "zero" on 
the spatial scale) with the most luminous row of the CCD 
("central row"); this does not have to coincide with the 
dynamical center. NGC 3115 or NGC 4179 (both classified 
as SO galaxies in the NASA/IPAC Extragalactic Database 
(NED)) were used as templates to subtract the stellar 
contribution from each single row. These template spec- 
tra turned out to be rather similar in our spectral ranges 
and were usually suitable to fit the absorption lines of the 
megamaser galaxies. Only for IC 1481, a special treatment 
was necessary that will be discussed in the next section. 

For line ratios presented here, the three central rows 
were averaged to enhance the S/N without loosing any 
spatial information. Hence, the line ratios refer to the cen- 
tral region of 2" x 2". Forbidden-line wavelengths were 
taken from lBowenI l( 1 960? ) . Heliocentric corrections as given 
in TableQwere added to the observed velocities. The sam- 
ple properties are given in Table |21 

3. Results 

The dereddened spectra of the innermost regions (2" x 2") 
rebinned to rest wavelengths are shown in Fig. ^ 

3.1. Stellar population of the galactic nuclei 

To derive emission-line fluxes, we need to model the stel- 
lar population which can contaminate our emission-line 
spectra due to underlying absorption lines. As a flrst ap- 
proximation, we simply subtracted one of the SO-type 
template spectra scaled to the inter-emission-line parts 
from all spectra. For two of the three objects, there re- 
mains no notable continuum nor absorption-line spec- 
trum. This means that, to our accuracy, the stellar popu- 
lations in the nuclei of ESO 103-G035 and TXS 2226-184 
closely resemble the old stellar populations in NGC 3115 
an d NGC 4179. 

I.Toguet et all l|200lh classifled ESO 103-G035 as 
"Seyfert 2b" meaning that they observed high order 
Balmer lines in absorption which are often interpreted 
to indicate a signiflcant young stellar population. We 
do not find high order Balmer absorption lines in our 



^ Munich Image Data Analysis System, trade-mark of the 
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Table 1. Summary of the spectroscopic observations 



Object/hel. corr." 


Date (beg.) 


p. a. (slit) 


A range (A) 


Exp. time (s) 


ESO 103-G035 


14-May-96 


44° 


5400-7400 


1800 


+ 14km/s 


17-May-96 


44° 


3400-5400 


1800 


TXS 2226-184 


14-May-96 


50? 3 


5400-7400 


2700 


+29 km/s 


17-May-96 


50? 3 


3400-5400 


2700 


IC 1481 


14-May-96 


24° 


5400-7400 


3600 


+25 km/s 


17-May-96 


24° 


3400-5400 


2400 



" This heliocentric correction was added to the measured radial velocities. 



Table 2. H20-megamaser sample properties 





ESO 103-G035 


TXS 2226-184 


IC 1481 


literature 


coordinates" 


18''38""20= -65°25'39" 


22''29™13' -18°10'47" 


23''19""25= +05°54'22" 


incl.''(°) 


67° 


70° 


62° 


m.a.'^(°) 


44° 


51° 


70° (30°"^) 


Vopt" (km s"-^) 


3983 ± 34 




6118^ ± 50 


Type« 


Syl.9/2, SaO 


?, SO? 


LINER, Sb pec' 


this work 


Whci" (km s ^) 


4040 ± 20 


7520 ± 30 


6120 ± 30 


dist.' (Mpc) 


61.9 


110.6 


88.5 


lin. sc. (pc/") 


300.2 


536.3 


428.8 


(IO^Mq) 


1.2 






Type 


Sy2 


LINER 


LINER 



" J2000, NED 

' Inclination angle from face-on (RC3) 

Position an gle of major axis (R C3) 

Taken from lBraatz et alJ lll997ll 
" Heliocentric velocity from RC3 
f Taken from I Van den Bergh et"ai] (|20o2) 
^ Classification and Hubble type (NED) 

Heliocentric velocity as derived from our rotation curves. 
' Distance via direct application oi v = Ho ■ d {Ho — 65kms~ Mpc~ ) using velocities relative to the 3K background as 
derived from Whoi with the velocity calculator provided b y NED . 

Estimated lower limit of mass with q=0.6 (.Lcaucux> Il983r within the central 2.8 kpc (diameter). 



spectra. Thi s migh t be due to our lower S/N spectra: 
iJoeuet et a D ('2001') used the same telescope and instru- 
ment but a UV -blazed grating and a CCD camera with 
high quantum efficiency below 5000 A. Furthermore, they 
used twice the integration time with half the spectral reso- 
lution. However, having a closer look at their online data, 
we find that E SQ 103-G035 resemb les more closely some 
galaxies which I.Tog^iet et al ] (|200lh classified as "Seyfert 
2a" (spectra showing both emission and absorption lines 
in which the Balmer series is not detected in absorption, 
e.g. ESO 323+G32) than "Seyfert 2b" , thus in agreement 
with our spectra. In our raw data we see Can KA3933, 
weak CNA4200, G bandA4301 and Mgi A5175 absorption 
lines. Since all of them can be well fit by NGC3115, we 
exclude the possibility of a significant residual starburst 
contribution. 

For IC1481, the subtraction of SO-type spectra was 
not satisfying as strong Balmer absorption lines and weak 
Can H+K features remained. Thus, we modelled the un- 



derlying stellar continuum in two ways: i.) by using the 
PEGASE code (version 2.0) bv lFioc fc Rocca-Volmerangel 
(1997) and ii.) by using the best fit of linear combinations 
of main sequence A to F stars and the template spectrum 
of NGC 3115. The latter provided a much better fit to the 
observed spectra of IC1481, and was thus used. The fit 
was performed as follows: 

We tried a variety of linear combinations of main se- 
quence A t o F star s on t he one hand (taken from the 
catalogue of lPicklesI l|l998fl . rebinned with a spline func- 
tion to the higher resolution of our spectra), and the tem- 
plate spectrum of NGC 3115 on the other hand to fit 
the blue spectrum with its strong features and simulta- 
neously the spectral distribution in the red range. The 
scaled monochromatic fluxes at 5556 A were weighted ac- 
cording to the different fluxes of AOV — F8V stars, and, 
by using a Salpeter initial mass function, normalized to 
the flux and mass of an A5V star. 
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ESO 103-G35 



i 300 



I 




x[A] 

TXS 2226-184 




|(l)yi|^^f^%y.#w#y«^^'' 



5000 6000 
X[k] 



Fig. 1. Combined dereddened spectra of the red and 
blue wavelength range rebinned to rest wavelengths for 
comparison. A template galaxy spectrum (NGC 3115 
or NGC 4179) was subtracted for ESO 103~G035 and 
TXS 2226-184. The upper observed spectrum of IC 1481 
shows strong Balmer absorption lines and weak Call H+K 
features. Thus, a linear combination of AOV — F8V stars 
and the template galaxy was subtracted (see middle spec- 
trum). The result is shown in the lower spectrum. (Both 
upper spectra of IC 1481 are shifted vertically by an arbi- 
trary amount.) 

After a first approximation of such a linear combina- 
tion to gain the absorption-line free Ha/H/3 value, the 



spectra of IC 1481 were dereddened using the recombina- 
tion value for the intensity ratio Ha/H/J — 3.1 (a typ- 
ical value for AGN s) and an average reddening curve 
ijOsterbrockl l(l98!^ . Table 7.2; MIDAS command "ex- 
tin/long"). This procedure was carried out several times 
to find the best combination iteratively (E'b-v = 0.15): 



fic 1481 = 0.099 /aov + 0.076 /a2V + 0.076 /a3V 
+0.071 /a5V + 0.064 /a7v + 0.061 /fov 
+0.039 /f5V + 0.036 fpsv + 0.48 /ngc3ii5 



10 ^^ergs ^ cm 



(1) 



As the lowest spectrum in Fig.^shows, there are no signif- 
icant absorption features left that could constrain a more 
elaborate population model to be subtracted. A slight de- 
pression near ^ 4020 A might be suggestive of Hel A4026 
from early B stars but is too noisy for detailed fitting. 
There is no evidence for He I A4471 so that we estimate a 
main-sequence turnoff at late B stars. This yields a rough 
age of an a s sume d instantaneous starburst of ^ 1.7 • 10* 
years lllbenl . ll967l) . 

We converted the theoretical physical continuum fluxes 
of the above stars as tabulated by Cox ( 2000) into fluxes 
which would be measured at the assumed distance of r = 
88.5 Mpc (IC1481, Table EJ using fx = Fx^ {R = 1.2 
— 2.9 Rq) and compared them with the measured ones. 
This leads to 9.4 -10^ stars or a mass of 1.7 • 10* Mq in the 
range AOV to F8V. In a more complete range for the orig- 
inal starburst (05V — M5V, mass range 60 — 0.21 Mq), 
this would correspond to a total mass of newborn stars of 
9.9 • 10* Mq or 1.4 -10^ stars. Although these estimates 
are rather crude, they nevertheless show that there was a 
vigorous starburst in the central kiloparsec of IC 1481 10* 
years ago. This finding is contradictory to iBoisson et alJ 
|2000) who claimed that in their sample of central re- 
gions of 12 galaxies (normal galaxies, starburst galaxies, 
LINERs, Seyfert Is & 2s), LINERs show the oldest popu- 
lation with little star formation still going on. It is, how- 
ever, not inconsistent with a n HST image that was an- 
alyzed bv iFalcke et all l|200l|) . They speculated that the 
irregular shape of the galaxy might be the site of an on- 
going galaxy merger. 

It could be interesting to perform a detailed empir- 
ical popula ti on synthesis as propo s ed and developed by 
iBical iTigsi. ISchmitt et a,lJ lT99i. ICid Ferna,ndes et all 
l|200ll) . and ICid Fernandes et alJ l|2003l) for our three 
galaxies and to compare the results with other stellar 
populations o f Seyfe r t-2 and LINER-host galaxies (e.g. 
ISchmitt et all l|l999t) : iBoisson et all l|200nt) : I.Toguet et alJ 
f2001)). However, on the one hand, having a very small 
sample (one Sy2 and two LINERs), we could not provide 
deeper insights into the links between nuclear activity and 
the star formation history of the host galaxy. On the other 
hand, as we do not expect to have missed a significant con- 
tribution of underlying stellar absorption lines, this time- 
consuming procedure is beyond the scope of this work. 
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3.2. Rotation curves 

To gain velocity curves, we took averaged peak wave- 
lengths of Ha and [N il] A6583 (with heliocentric correc- 
tion applied, see Tabled) and subtracted the hehocentric 
systemic velocity derived by symmetrising the curve of 
each galaxy. This heliocentric velocity can be found in 
Table [3 Rotational velocities were calculated assuming 
that all spectra were taken along the major axis of each 
galaxy. Thus, the observed velocities had to be corrected 
simply for the inclination angles (as taken from RC3): 



Rotation curve of ESQ 103-G35 



Vrot 



Vphs 

sini 



(2) 



For the central 2-4 kpc for which the Ha and [N ii] A6583 
lines were strong enough to be measured, the rotational 
velocities are presented in Fig. [3 The central mass could 
only be estimated for ESO 103~G035 (see Table 13, be- 
cause the "rotation" curves for IC 1481 and TXS 2226-184 
are remarkably flat. As a lower limit, we estimate the mass 
as 



M = 



■ R 



G 



(3) 



with a = 0.6, taking into account non-Keplerian motion at 
those locations, where the gravitational potential consists 
of the superposition of a flat disk and a spherical com- 
ponent (LcQucux, 1983). The resulting mass of ESO 103- 
G035 is 1.2 X 10^ Mq within the central 2.8 kpc {Tahle^. 

The flat "rotation" curve of IC 1481 can be explained 
by its amorphous appearance: The spectra were possibly 
not observed along the major axis, since the position of 
the major axis is very uncertain. The observations were 
done at a position angle of ^ 24° correspo nding approx- 
imatel y to the position angle reported by iBraatz et alJ 
However, in the RC3, the position angle of the 
major axis is listed as 70° (Table [^J. 

TXS 2226-184 was definetely observed along the major 
axis. This galaxy might therefore be an elliptical with its 
central gas disk lying coplanar to the major axis. While the 
latter is strenghtened by observations of the narrow-line 
region, the surface brightness distribution favors a spiral 
galaxy ( Falckc ct al... 2000'). The flat rotation curve may 
thus be caused by the fact that our observations arc con- 
fined to the extent of the nuclear bulge where gas kinemat- 
ics do not follow the trend expected for a highly inclined 
disk that should be dominant further out. 



north-east Distance [pc] south-west 

Rotation curve of TXS 2226-184 



north-east 



Distance [pc] 



Rotation curve of IC 1481 



south— west 



Distance [pc] 



south-west 



Fig. 2. Rotation curves of the three H2 0-megamaser 
galaxies derived from the average value of the Ha and 
[N 11] A6583 peak wavelengths. 



3.3. Emission-line profiles, intensities and widths 

Attempts to flt the emission- line profiles of the strongest 
fines from H, [OIII], [01], [NH], and [SII] in the tem- 
plate subtracted spectra were made with Gauss functions, 
Lorentz functions (— Cauchy distributions), and modified 
Lorentzians, dubbed here as d-Lorentzians, as described 
in Appendix 1X1 (Note, that in the MIDAS environment, 
Lorentz functions are called Cauchy functions, while d- 
Lorentzians are called Lorentz functions.) 



It turns out that all line profiles from the galactic nu- 
clei had stronger wings than a single Gaussian and slightly 
weaker wings than a single Lorentzian, but could usually 
be well fit by a d-Lorentzian, which also allowed an appro- 
priate decomposition of the Ha+[NII] blend. In Table 13 
heliocentric velocities and line widths (FWHM) as well 
as mean d-parameters as derived from d-Lorentz fits are 
shown for the six strongest lines (excluding [0 11] A3727, 
as it is an unresolved double line). Examples of the fits 
are shown in Figs. 01 and ^ 
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Table 3. Heliocentric velocities and line widths (FWHM) (separated by |) in kms ^ as well as mean d-parameter as 
derived from d-Lorentz fits for the central row. 



line 


ESO 103-G035 


TXS 2226-184 


IC 1481 


H/3 


4011 1 479 


7498 


446 


6138 1 445 


[0 III] A5007 


4045 539 


7553 


397 


6123 405 


[O i] A6300 


4063 584 


7564 


676 


6087 405 


Ha 


4046 504 


7486 


327 


6098 392 


[Nil] A6583 


4046 443 


7486 


526 


6097 416 


[S II] A6731 


4059 1 411 


7515 


389 


6107 1 363" 


mean 


4045±8|493±26 


7517±14 


460±51 


6108±7 1 404±11 


mean d 


1.6±0.2 


l.lzbO.l 


1.4±0.3 



" [S ii] A6731 was truncated by telluric absorption bands and thus, [S ii] A6716 was used instead. 
Table 4. Observed and reddening-corrected line intensity ratios relative to H/3.° 



Line ESO 103-G035 TXS 2226-184 IC 1481 





Fobs 


-fdcrcd 


-fobs 


-fdcrcd 


-fobs 


-fdcrcd 


[0 II] A3727 


1.53 


2.56 


4.39 


5.69 


3.03 


3.43 


[Ne III] A3868 


0.67 


1.03 


0.4 


0.49 


0.51 


0.57 


He+ [Nelll] A3967 


0.22 


0.33 










[0 III] A4363 


0.12 


0.15 










HciiA4686 


0.14 


0.15 


0.18 


0.19 






[0 III] A5007 


9.93 


9.44 


2.81 


2.74 


1.04 


1.02 


[Fell] + [Fevil] A5159 


0.13 


0.12 










[Fevi] A5176 


0.05 


0.04 


0.14 


0.13 


0.13 


0.12 


[Ni] A5199 


0.18 


0.15 


0.3 


0.28 


0.31 


0.3 


HeiA5875 


0.28 


0.18 


0.58 


0.47 






[0 i] A6300 


1.62 


0.93 


2 


1.42 


0.96 


0.83 


Ha 


6.11 


3.1 


4.45 


3.1 


3.66 


3.1 


[Nil] A6583 


6.4 


3.25 


5.13 


3.57 


2.36 


2 


[Sii] A6716 


2.14 


1.05 


1.56 


1.05 


1.5 


1.26 


[S II] A6731 


2.1 


1.03 


1.48 


0.98 






HeiA7065 


0.06 


0.03 










[Arm] A7136 


0.49 


0.22 











" Derived from the central region (2" x 2") 



Table 5. Reddening-corrected H/3 luminosity and results from dereddened line ratios. 





ESO 103-G035 


TXS 2226-184 


IC 1481 


iH,3" (10-*"ergs~^) 


2.65 


0.49 


0.59 


Av" 


1.9 


1. 


0.47 


ne" (cm-3) 


600 ± 50 


550 ± 100 




T*" (K) 


13000 ± 500 






[Oiii] A5007/[Oii] A3727 


3.69 


0.48 


0.3 


[0 III] A5007/[O i] A6300 


10.15 


2.08 


1.23 


[Oi] A6300/Ha 


0.3 


0.46 


0.27 


[Sii] A6716-f A6731/Ha 


0.67 


0.66 


0.79" 


[Nil] A6583/Ha 


1.05 


1.15 


0.65 



" Derived from the central region (2" x 2") 
' Derived from the central row (0'.'68 x 2") 

° Using an estimated ratio of [Sll] A6716/A6731 ~1.07 for a typical ric of ~ 550 cm 

From the parameters of the d-Lorentz fits of strong 
lines, the line intensities relative to H/3 were derived (see 
Appendix ^ while the intensities of faint (single) lines 
were directly obtained by integration using the MIDAS 
command "integrate/line". As for IC 1481, extinction was 



derived with the c-method IjOsterbrockl Il983) by using 
the recombination value for the intensity ratio Ha/H/3 = 
3.1 and an average reddening curve. Both observed and 
reddening-corrected line-intensity ratios from the central 
region (2" x 2") are presented in Table0] For pairs of lines 
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([O III], [O i], and [Nil]) with a fixed fine ratio (~ 3:1), only 
the brighter line is fisted. Estimated errors of the line ra- 
tios lie in the range of 10% — 20%. Comparison with liter- 
ature was only possible for th e brighter lines of ESO 103- 
G035 ijMorris fc Wardl.ll988(l . as all other fine ratios have 
not been reported in the literature yet. Deviations of ob- 
served line ratios lie in the ran ge of ^ 1 % — 35 % and 
may be due to the fact that Mo rris fc WardI l|l988|) did 
not use a template spectrum to subtract the underlying 
stellar contribution. Furthermore, they refer to a larger 
region (1'.'5 x 5"). 

Applyin g the classical methods outlined in 
lOsterbroc^ l|l989l) . we derived the electron density 
using [Sli] A6716/A6731 in the central region (2" x 2"). 
The temperature, determined from [Oiii] ^^^^^§^, 
only corresponds to the brightest row (formally O'.'BS x 2") 
because of the weakness of the A4363 fine. Clearly, these 
parameters given in Table represent some average 
over the central several hundred parsecs in each case 
(see Table |3 for the typical scale in each galaxy). For 
TXS 2226-184 and IC1481, no temperature was mea- 
sured, since the [Oiii] A4363 emission line was too faint. 
In the spectrum of IC1481, the [S ii] A6731-emission 
line was truncated by telluric absorption bands, thus no 
electron density could be derived. For TXS 2226-184, 
[S ll] A6716 might be slighly affected by the telluric 
absorption bands and the deduced electron density has 
to be taken with some caution. The reddening-corrected 
H/3 luminosity for the central region (2" x 2") is given in 
Table El as well. 

The most important fine ratios from the central re- 
gion (2" X 2") which distinguish between emission-line 
object classes were calculated (se e Table and pu t in 
the diagnosti c diagrams as used in lOsterbrockl l|l989tl and 
iBaldwin et a l. (1981). ESO 103-G035 falls into the Seyfert 
regime. As the spectra of ESO 103-G035 did not show 
any broad components neither in Ha iior in H /3 in agree- 
ment with the results of lMorris fc WardI ljl988f) . we classify 
ESO 103-G035 as Seyfert-2 galaxy. 

Line ratios of TXS 2226- -184 and IC 1481 follow closely 
the LINER classification of lHeckmanI 



Gauss Fit to Hf, [OIII] (ESO 103-035) 



([Oiii]/[Oii] 

< 1 and [Oiii]/[Oi] < 3) and fall into the LINER regime 
in b oth the line-r atio diagrams used bv lOsterbrockl l)l989l) 
and lBaldwin et al . (1981). Furthermore, this classification 
is in agreement witfi their low l uminosity, and , for IC 1481, 
identical with the results from lHuchra et till l(l982f) . Our 
classification of TXS 2226-184 as a LINER is new. 



4. Discussion 




MA] 

d-Lorentz Fit to H/S, [OIII] (ESO I03-G35) 




Fig. 3. H/3 and [O ill] AA4959,5007 fit by three Gaussians 
(upper), three Lorentzians (middle) and three d- 
Lorentzians (lower) (ESO 103-G035). The total fit is 
shown (dashed line). Gaussians yield too narrow wings, 
Lorentzians too broad ones, whereas d-Lorentzians give 
the best total fit. 



4.1. Why d-Lorentzians? 

When discussing emission-line profiles of active galaxies, 
it is important to keep in mind, that they might represent 
line-of-sigfit integrations of several kinematical compo- 
nents such as re stating disks, (jOnes of outflowing gas or ex- 
panding shells. 'Schul z et all ( 199!tIi have shown that line- 
of-sight integrations through such configurations would 



neither yield Gaussians nor Lorentzians. Furthermore, if 
the spatial and spectral smearing by an instrumental func- 
tion is wide compared to the kinematical gradients, the 
profiles will reveal a more symmetric core-wing structure 
l|Schulz et allll995l) . The latter also applies here with our 
spectral resolution of ~ 130kms~^ and spatial resolution 
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Gauss Fit to Ha. [Nil] (IC 1481) 




«: 10 



some kind of turbulence, they often fail to fit broad wings. 
They can instead be well fit by Lorentzians, which dif- 
fer fro m Gaussians w ith same FWHM by more extended 



Lorentz Fit to Ha, [Nil] (IC 1481) 




d-Lorentz Fit to Ha. [Nil] (IC 1481) 




A[S] 

Fig. 4. The same as in Fig. |3 for the Ha and 

[Nil] AA6548,6583 blend (IC1481). Again, d-Lorentzians 
lead to the best total fit. 



of ~ 1'.'2 - 2'.'5 which lead to an integration over relatively 
large cells of the phase space of an NLR. 

Various observations revealed that even the "narrow" 
lines in AGNs can show considera ble profile structure 
if measured at sufficient resolution ( Vrtilek fc Carletonl 
ll985HWhitti^ . ll985tlsd^ulz fc Hen kell . 1200.'^ . Commonlv. 
in the case where single-component fits fail, multi- 
component Gauss fits or Lorentz fits are employed. 

While Gaussians are relatively well suited to describe 
instrumental functions, thermal Doppler broadening or 



winefs. lwhittlel ^^198{^^ already described the non-Gaussian 
nature of observed [O III] line profiles which revealed 
"a stronger base relative to the core than Gaussians" . 
IVeron-Cettv et aP l|200lj) also reported that most broad 
emission lines of Narrow-Line Seyfert 1 (NLSl) galaxies 
could be well fit by a single Lorentz profile, confirming 
previous claims that Lorentzians rather than Gauss pro- 
files are better suited to reproduce the s hape of the NLSl's 
broad emission lines. ISchulz fc Henkei l|2003l) showed that 
Lorentzians appear to be better suited to fit the spectra 
of NGC 1052 and Mrk 1210. Hence, extended wings of in- 
trinsic narrow-line profiles of the bulk components may 
not be uncommon. 

However, the use of d-Lorentzians is new and allows us 
to fit both permitted and forbidden lines by adjusting the 
additional parameter d (see Appendix for consequences 
concerning the FWHM). Note that, as we are dealing here 
with Seyfert-2s and LINERs, no broad lines were fit. In 
those cases in which both broad- and narrow-line com- 
ponents can be clearly distinguished, a multi-component 
fit is of course justified. According to Occam's razor, the 
use of d-Lorentzians has the advantage of adding just a 
single parameter instead of using multi-component fits of 
Lorentzians and Gaussians. As d-Lorentzians seem to pro- 
vide a very good fit to the data, the emission-line fiuxes 
can be well approximated as the integral over the corre- 
sponding d-Lorentzian, which is a great advantage espe- 
cially when measuring blended lines. Using d-Lorentzians 
does not imply a loss of physical information: On the one 
hand, both Lorentzians and Gaussians might not be phys- 
ically motivated in all cases (as argued above) and on the 
other hand, the parameter d gives a measure of the width 
of the line wings. 

Our observed line wings have typically reached the 
noise levels at 6% of the maximum. Using equations 
(A.l) — (A. 3), a Gaussian falls to 6.3% of the maximum at 
2 HWHMs (half width at half maximum) from the center, 
a Lorentzian to 5.9% at 4 HWHMs and a d-Lorentzian has 
decreased to 5.9% at 2.8 HWHM, taking a typical value of 
d ~ 1.3. Using a mean FWHM from the six strongest lines 
of the three galaxies (Table O)), the line wings have veloc- 
ities of 600 — 700kms~^. Thus, the question arises how 
such large velocit ies in the line wings can b e attained. This 
was discussed by ISchulz Hen kdl ll200,'^ in terms of the 
presence of turbulence, outflows, magneto-hydrodynamic 
waves and electron scattering. They concluded that the 
latter might provide a viable explanation of the observed 
line wings. 

For our new sample of galaxies, we may ask again, 
whether electron scattering can produce the observed line 
wing velocities, cither due to the presence of hot coronal 
or cooler ionized gas. The ratio of scattered to input lu- 
minosity can be estimated as Ls/Lin = /te, with Tg = 
electron-scattering optical depth and / = covering frac- 
tion of the scattering medium. To obtain Te on the per- 
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cent level, effective column-densities of fully ionized gas of 
2-10^^ cm~^ would be required. The total bremsstrahlung 
luminosity of a kpc-extended sphere of a hot intercloud 
medium with such a column density can be estimated as 
2.6 ■ 10^ ^(n./10cm-^)\/r/10 6K(i?,/500pc)3erKS-i (Eq. 
5.15b in iRvbicki fc LightmanI |1979D). The intrinsic (un- 
absorbed) thermal soft -X-ray emission for ESO 103~G035 
(the only galaxy in our sample with measu red 0.5-2 keV 
flux) is estimated to be 4 • lO^^ergs"^ 



I Turner et all 



^23) ■ For a temperature ofr~ 10^, a density of ne ~ 10^ 
and a radius of i? ~ 1 kpc, a lumi nosity ^ 500 times higher 
would be expected. Thus, as in ISchulz & He^ ^l2003^ 
for NGC 1052 and Mrk 1210, we do not see enough hot 
coronal gas to produce the required wings in ESO 103- 
G035. However, electron scattering inside the cooler ion- 
ized gas component (T ~ 10^ K) with electron densities 
of Ue 500 (which were observed for ESO 103-G035, 
Table 01 is a possible explanation. Calculating the ex- 
pected soft-X-ray luminosity with a scale of i? ~ 10^ pc 
yields ~ 5- lO^^ergs" ^ . in a greement with the measure- 
ments of iTurner et al.l l|l997t) . Considering a geometrical 
scattering factor / ~ 0.3, one easily obtains Ls/Lin of a 
few percent, which suffices to explain the extended wings 
of the narrow lines. Thus, electron scattering even at the 
cooler ionized gas can account for the observed velocities 
in the line wings. 



4.2. What triggers the water megamaser activity? 

The unique association of H2O megamasers with 
AGNs of Seyfert-2 and LINER type and the fact 



sec(s) 


of the 


parent 


galaxy (IGreenhill et al.. 


1995a 
et al. 


Mivoshi ct al.', 'l995'; 


'Clausscn et al. , Il998l iTrottei 


199& .Herrnstein et al.,. 


.1999,) suggest that the unknown 



excitation mechanism is closely related to the AGN ac- 
tivity (and that at least some LINERs are AGNs rather 
than starbursts). The Seyfert-2 geometry is favorable for 
megamaser activity supplying the necessary large column 
densities of warm and dense molecular gas enriched with 
H2O molecules. 

Not in all cases, however, a nuclear radio contin- 
uum background has to be a mplified to obtai n detectable 
megamaser emission (e.g. IGreenhill et all lll995albl) 'l. 
Sometimes, column densities of a disk may be large enough 
even when viewed face-on and little is known about the 
solid angle of the emitted maser radiation. It is thus pos- 
sible that in exceptional cases nuclea r maser emis s ion is 
also detectable in Sey fert-1 galaxies l|Nagar et all l20o3 
iHagiwara erall 12003(1 . 

Interaction of molecular gas with radio jets will be 
associated with shocks an d the shock heate d gas favors 
II2O maser emission (e.g. IPeck et all l|2003|) 'l. Optically 
detectable ionized o utflowing gas may, h owever, also trig- 
ger maser emission. iFalcke et al.l l|200l|) reported linear 
jet-like HST features in TXS 2226-184 and NGC 1386 
that might be related to an outflow rather than an exci- 



tation cone. In the megamaser sample of lSchulz fc Henkel 
(2003), galactic rotation and outflow of narrow-line gas 
are common features. The detection of H2O megamasers 
in very luminous infr a red galaxies (Lf tr > 10^^ L©; 
IHagiwara et~al] l|2002al) : IPeck et alJ l|2004^ ) with presum- 
ably many young clusters of O and B stars may add fur- 
ther support for such an outflow scenario. 

4.2.1. Three specific cases 

There are a few nearby H2O megamaser sources that are 
particularly well studied, not only in view of their H2O 
maser properties but also at other wavelengths, including 
the optical and near infrared windows. Best known are the 
three northern of the five originally dis covered rtiegam aser 
sources, N GC 1068, NGC 4258 fciaussen et al.Lll984l) an d 
NGC 3079 llHenkel et al.L Il984l: iHaschick fc Baani Il985l) . 
The H2O data from NGC 4258 show a warped nu- 
clear accretio n disk that is seen almost edge-on (e.g. 
iMivoshi et all l|l995 ^). Towards NGC 1068, we may view 
a torus with "sub-Keplerian" rotation curve l|HureL l2002|l 
that might be more massive, also with respect to the nu- 
clear engine, and that appears to be thicker or slightly 
less well ordered. To the north, where the nuclear jet 
is ben ding, jet-masers are observed (e.g. ICalhmore et alJ 
1(200 ill ). NGC 3079 appears to show a disk that is less 
wefl ordered than those of NGC 1068 and NGC 4258, but 
a detailed map of the red-shifted m aser features i s still 
needed to obtai n a complete picture ((Trotter et all liggst 
Hagiwara ~t all |2D02b) . 

Direct hints for the presence of megamaser emission 
in these sources is given by the near infrared contin- 
uum, providing evidence for the presence of warm dust 
heated by the AGN. NGC 3079 harbors a compact (1") 
nuclear source wit h dust at a temperature of ~1000K 
((Israel et al.L Il998ll. An opticall v detected large scale out- 
flowing zone ((Cecil et al.l 1200 1() is apparently not trigger- 
ing maser emission. An even more compact dusty core 
is detected in NGC 4258, with a n upper size limit of 
200 mas (7pc; ICharv et aP ((200(tl ). Two jets (anoma- 
lous arms) reach out of the nucle ar region and shoc k 
the ambient gas in the inner 350pc ((Wilson et al.U200l() . 
Interestingly, the jets are much more prominent than in 
most other spirals, but do not trigger any megamaser 
emission. Instead, they provide the radio continuum back- 
ground for the enhanced flux of the systemic H2O features 
that are stronger than the presuma bly self-amplified red - 
and blue-shifted components (.Herrnstein et all Il998l). 
The b est studied galaxy is NGC 1068 re.g. lGalliano et alJ 
((2003!) ). This prototypical Seyfert-2 galaxy also contains 
a compact near- and m id-infrared source I Marco et all 
h997HRoua,n et all 119981) . The extremely red colors of its 
200 mas core (15 pc) lead to an intrinsic extinction of Ay 
25™, assuming classical grains at 1500 K. Optical spec- 
troscopy with the HST shows lines split into two veloc- 
ity compon ents separat e d by 1500 km s~^ within the in- 
ner arcsec ((Axon et all ^2^. 2 micron H2 emission is 
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also double peaked l)Galliano fc AlloiiJ. lioo^ . Both can 
be explained by gas interacting with the radio jet. 

To summarize, the three galaxies NGC 1068, 
NGC4258, and NGC 3079 exhibit a spatially com- 
pact near infrared core containing dust clouds that are 
heated by the central engine. This appears to be the main 
hint for the potential presence of accretion disk masers, 
but these cores may be more difficult to detect than the 
masers. Tracers for jet masers are split lines at optical and 
near infrared wavelengths while outflow masers should 
show broad wings in lines arisin g from ionized gas as e.g . 
seen in NGC 1052 and Mrk 1210 (Schul z fc Henkej 120031. 
The nature of these outflows has, however, to be clearly 
identified to distinguish them from gas entrained by the 
much faster nuclear jets. 

4.2.2. Our data 

The new data presented here support the connection 
between H2O megamaser emission and the presence of 
Seyfert-2 and LINER nuclei. They do not, however, 
further elucidate the relationship between outflows and 
megamasers as they do not show signs of outflow in their 
optical spectra. In any case, establishing a connection be- 
tween optical and radio data is difficult keeping in mind 
the different angular scales involved (i.e. a few milliarcsec- 
onds for 1.3 cm interferometric maps and a few arcseconds 
for our study), so that optical and near infrared data taken 
with subarcsec resolution would be highly desirable. 

4.2.3. H2O maser classification 

We conclude that all 22 GHz H2O masers detected so far 
might be subdivided into five categories with rising but of- 
ten overlapping ranges of luminosity: (i) Masers from late- 
type stars, (ii) masers from star forming regions, (iii) nu- 
clear "outflow" masers in which outflows impinging onto 
dense molecular clouds may provide a suitable trigger for 
kilo- or megamaser emission (possibly occurring in star- 
burst galaxies, LINERs dominated by starburst instead of 
AGN activity, and type-1 AGNs) , (iv) "jet" masers where 
a direct interaction between the nuclear jet and the inter- 
stellar medium triggers maser activity and (v) "accretion 
disk" masers where the masers are aligned in sometimes 
warped disks or thick tori. 

5. Conclusions 

We analyzed optical spectra of the megamaser galaxies 
ESQ 103-G035, TXS 2226-184, and IC 1481 and presented 
rotation curves, line ratios, electron densities, tempera- 
tures, and Lufj. The successful line-profile decompositions 
of single lines as well as blends rest on d-Lorentzians with 
an additional parameter d to adjust the wings, rather 
than Gaussians or Lorcntzians as basic functions. Using 
d-Lorentzians, we can fit the extended lines of intrinsi- 
cally narrow-line profiles. Crude estimates suggest that 
electron scattering at the ionized gas itself might lead to 



a viable explanation of the large line widths observed, 
but detailed modelling of such processes would be use- 
ful. On the basis of line ratios and their low luminosity, 
TXS 2226-184 and IC 1481 are classified as LINERs while 
ESO 103-G035 is a Scyfcrt 2. No significant velocity gra- 
dient is observed along the major axis within the inner 
2kpc of TXS 2226-184. IC 1481 reveals a post-starburst 
spectrum which could not be fit by an SO template galaxy 
only. Instead, using additional AOV — F8V star spectra, 
we could correct for the Balmer absorption lines. 

The comparison of optical spectra of all three 
H2O megamaser galaxies with an angular resolution of 
~ 2" support the connection between II2O megamaser 
emission and the presence of Seyfert-2 and LINER nuclei. 
However, the spectra do not provide direct clues to the 
very nuclear regions where the masers reside. Their clas- 
sification as "outflow", "jet" or "accretion disk" masers 
remain uncertain. Apparently, either the nuclei are too 
heavily obscured by dust associated with cool dense fore- 
ground gas or the morphological structures giving rise to 
the megamaser emission, i.e. accretion disks and jets, have 
linear scales that are far too small to be detectable with 
an arcsec resolution. 
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Appendix A: Line-profile fits 

Line profiles of the strongest lines have been fit with 
Gaussians, Lorentzians, and d-Lorentz functions as pro- 
vided by MIDAS using the same initial parameters and 
the commands "edi/fit" and "fit/image". These basic fit 
functions are compared in Fig. lA.ll where it is shown that 
Lorentzians have much stronger wings than Gaussians, 
but that the wings of d-Lorentzians can be adjusted to 
intermediate wing strengths. The functions are defined 
as (a — height, b — position, c = FWHM for Gaussian 
and Lorentzian, and d = additional parameter for a d- 
Lorentzian) 



Gauss {x; a, b,c) — a exp 



Lorentz (x; a, 6, c) 



d— Lorentz (x; a, 6, c, d) 



-In 2 



2{x-b) 



2ix - b) 



2{x-b) 



(A.l) 



(A.2) 



(A.3) 



Bennert et al.: Spectral characteristics of three water megamaser galaxies 



11 



In all galaxies, all spectral rows and for all analyzed lowing integral 

-^(i— Lorcntz — CL • C • 



Gauss-, Lorenlz-, and d—Lorentz— Function 




r(-o.5 + rf) 

2 f(d) 



(A.6) 
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Fig. A.l. Upper panel: Comparison of Gaussian (solid), 
Lorentzian (dotted) and d-Lorentz functions (dashed line) 
(a = 100, b = 5000, c = 20 for Gaussian and Lorentzian, c 
= 26.1 and d = 1.5 for d-Lorentzian to gain the same effec- 
tive FWHM). A d-Lorentzian falls in between a Gaussian 
and a Lorentz function. Lower panel: Four d-Lorentzians 
are shown with different d parameters [d — 0.8 (dashed), d 
= 1.2 (dotted), d = 1.5 (solid) and d = 2 (dash-dotted)]. 
With increasing d, the wing profile narrows. Thus, the pa- 
rameter d can be used to adjust the wings. 

strong emission-lines, a Gaussian yields too narrow wings, 
a Lorentzian too broad ones and a d-Lorentzian leads to 
the best fit (for comparison, two examples are plotted in 
Figs. 131 and 01). In a d-Lorentzian, c is not the FWHM 
width, which is instead given as a combination of c and d. 



FWHMd-Lorentz = C • . (A.4) 

The width at n % of the height can be calculated with 



Lorentz 



1/d 



(A.5) 



With increasing d the wing-profile narrows (see Fig. IA.1() . 
The line intensities can be derived by computing the fol- 
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